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. Abstract 
o\ . 

^ \ Transition from regular to chaotic dynamics in a crystal made of singular 

CIh' scatterers U{r) = Alrl"*^ can be reached by varying either a or A. We map 

' the problem to a localization problem, and find that in all space dimensions 

the transition occurs at o" = 1, i.e., Coulomb potential has marginal singu- 
larity. We study the critical line a = 1 by means of a renormalization group 
technique, and describe universality classes of this new transition. An RG 
' equation is written in the basis of states localized in momentum space. The 

^ ' RG flow evolves the distribution of coupling parameters to a universal station- 

ary distribution. Analytic properties of the RG equation are similar to that 
of Boltzmann kinetic equation: the RG dynamics has integrals of motion and 



o 



\ obeys an i?— theorem. The RG results for a = 1 are used to derive scaling 



laws for transport and to calculate critical exponents. 



V ■ I. INTRODUCTION 

Q \ In the theory of localization to which I. M. Lifshits made many seminal contributions [|l|, 

one is interested in the effect of random potential on the quantum-mechanical wavefunction. 
When the randomness is weak, the wavefunctions are extended throughout the whole system, 
^ . whereas at sufficiently high disorder all wavefunctions become localized 0. The transition 
between localized and extended states for many years has been at the center of interest in 
condensed matter physics . Originally, the concept of localization was invented in order 
to understand properties of disordered metals and doped semiconductors. More recently, it 
became clear that the ideas and the language developed in the localization theory are useful 
in a much broader context. One of the areas of theoretical physics that has numerous links 
with the localization theory is Quantum Chaos. In this article we consider quantum chaos 
in a periodic Kepler problem, and solve this problem by using tools from the localization 
theory. The relation between the Quantum Chaos and Localization theories appears to be 
quite generic, and we believe that it will help to understand the transition from integrable 
to chaotic behavior in a broad class of quantum systems. 

Traditionally, the field of Quantum Chaos deals with the relation between classical and 
quantum dynamics in chaotic systems. One is interested in the features of classical dynamics 
that are manifested in quantum dynamics, such as periodic orbits, ergodicity, etc. Another 
group of questions concerns the transition from integrable and chaotic quantum dynamics. 
Already in the classical case this transition appears to be a difficult problem, although 
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tractable within the KAM theory On the other hand, the transition in the quantum 

case, apparently, is not quite so well understood. For the reader's sake, below we review 
some basic facts. 

It is known that classical Hamiltonian dynamics can be either integrable or chaotic. 
Integrability means that all variables can be separated, which leads to quasiperiodic motion, 
concentrated in a small part of the phase space. In contrast, the chaotic (nonintegrable) 
dynamics is ergodic, i.e., the system can be found arbitrarily close to a given point in the 
phase space. 

One of the most challenging problems in the theory of classical dynamical systems is the 
transition from integrable to chaotic behavior when an integrable system is perturbed, and 
the perturbation violates the integrability. There exist two scenarios of this transition. In 
the case of a singular perturbation, even a perturbation of an arbitrarily small amplitude 
can drive the system into chaos. A smooth perturbation acts in a different way: according 
to the seminal KAM scenario (see the motion remains non-ergodic under a smooth 

perturbation if the latter is weak. 

A good example of a singular classical Hamiltonian dynamics is the motion in a peri- 
odic lattice of Coulomb scatterers 0. The dynamics in this problem can be mapped onto 
the problem of geodesies on a negative curvature compact space, which in the dynamical 
systems classification is known to be of the most unpredictable, and yet deterministic kind. 
The degree of ergodicity one has in this problem is expressed mathematically by so-called 
Bernoulli symbolic dynamics scheme. This dynamics leads to very fast mixing in the phase 
space quantified by Kolmogorov entropy. Another interesting example of a similar kind is 
the anisotropic Kepler problem |]^: a particle with anisotropic mass tensor moving in the 
Coulomb potential. It was shown by Gutzwiller that the existence of the trajectories that 
fall on the center breaks integrability at any potential strength. Qualitative features of the 
dynamics are similar to that of the periodic system of Coulomb scatterers, although the 
Bernoulli scheme is constructed in a different way [|^. 

In the quantum problem, a lot of interest recently was focused on the semiclassical 
theory of quantum billiards. The billirads represent one of the simplest problems in which 
one can study the relation between classical and quantum chaos (see review in The 
quantum problem of a periodic array of Coulomb scatterers was considered by Gutzwiller 
who developed a scattering mode representation well suited for treating the semiclassical 
limit. 

In this article we consider the integrability to chaos transition in quantum systems, focus- 
ing on the case of singular perturbation. The usual approach is to look for critical strength 
of the perturbation. Instead of that, however, we found it more natural and interesting to 
add another dimension to the problem, and to characterize the integrable to non-integrable 
transition as a function of the degree of singularity of the perturbation. It turns out that if 
the perturbation is both not strong and not singular enough the quantum behavior, namely 
the spectral statistics, will remain similar to the integrable one. On the other hand, when 
either the perturbation itself or its singularity is sufficiently strong, the behavior is chaotic. 
We also find an intermediate behavior for weak perturbations with the singularity taking 
certain critical value. This behavior can be called quasichaos, or weak chaos. 

The problem we will be interested in is the motion of a quantum particle in a three- 
dimensional crystal of scatterers with Coulomb cores. It will be assumed that the particle 
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velocity is so high that the wavelength is much smaller than the crystal period. This problem 
describes transmission of fast charged particles through a crystal, assuming the particles 
energy is so high that scattering by the atomic core potential is important. (We ignore 
complications arising from inelastic scattering.) The relation of our problem to the real 
situation corresponding to experiments on channeling particles in crystals will be discussed 
at the end of the paper in Sec. |X[ 

Usually, for a particle moving in a periodic potential there is a notion of a Bloch state, 
i.e., the particle wavefunction can be represented as a superposition of few plane waves. 
However, when the scattering by the cores is strong the number of plane waves forming each 
Bloch state becomes very large, and the wavevectors are distributed roughly uniformly over 
all directions in the momentum space. This means that the particle is not moving along a 
straight line, but rather is diffusing. These two regimes, ballistic motion roughly along a 
straight line and diffusive motion with an irregular change of direction, correspond to the 
Bloch states localized or delocalized in momentum space. 

The idea of our approach is the following. We treat the core potential as a perturba- 
tion, and write the problem using plane waves as a basis. In this basis, the main part of 
Hamiltonian, the kinetic energy, is diagonal. Since the scattering potential is periodic, the 
wavevectors of relevant plane waves form a lattice in momentum space. In the chosen basis 
we have a tight-binding localization problem: the potential matrix elements define hopping 
between different lattice sites, and the kinetic energy can be interpreted as on-site energy. 

For a non-singular core potential one effectively has only short-range hopping in the 
localization problem in momentum space. In this case, according to the conventional theory 
of localization, all states remain localized if the potential strength is below some threshold. 
On the other hand, if the core potential is singular, its Fourier components decay slowly 
and one has to deal with a localization problem with long-range hopping. We will see that 
long-range hopping canm delocalize the quantum states even at arbitrarily weak potential, 
if the potential singularity is sufficiently strong. 

The states localized in momentum space correspond to quasiperiodic motion in real space, 
i.e., to almost integrable dynamics, whereas the states delocalized in the momentum space 
correspond to chaotic dynamics in real space. Generally, the transition from one regime 
to the other can be reached in two different ways: either by increasing the strength of the 
core potential, or by varying the power law of the potential singularity at the origin. Below 
we consider the situation when the potential is weak: Ze'^/hv < 1, where v is the particle 
velocity, and Z is nuclear charge. In this case, naively, the scattering is weak and thus 
the Bloch states consist only of few plane waves. However, we will see below that because 
of the 1/r— singularity of the Coulomb potential the condition of weak scattering is more 
stringent: Ze'^/hvlnpaB/h < 1, where p is the particle momentum, and as is Bohr's radius 
of the cores. We will be interested in the situation when 

Ze^/hv < 1 , and Ze^ /hv \n{paB/h) > 1 , (1) 

In this case the scattering is strong, there is no ballistic motion, and the Bloch states consist 
of many plane waves. 

However, in this case the motion is not ergodic because the number of plane waves 
forming the Bloch state is much less than the total number of plane waves with appropriate 
energy. The reason is, as we will see below, that this problem falls right on the critical line 
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of a localization transition. To study this transition we introduce a more general model 
with core potentials having an arbitrary power law singularity (see Sec. |T|), and consider 
the transition from localized to extended states in momentum space resulting from changing 
the power law. The interesting thing about this transition is that it occurs even in the weak 
coupling limit, whereas in the conventional Anderson localization problem with short-range 
hopping all states are localized at weak hopping. As a result, the nature of the transition is 
quite different from the conventional picture accepted for short-range hopping, which leads 
to a completely different critical behavior. 

After introducing the model in Sec. |T|, and discussing the mapping to a localization 
problem in Sec. |ITI|, we proceed with solving the problem. It turns out that this problem is 
much more tractable than the conventional localization. The localization transition driven 



by long-range hopping has been studied previously and is well understood [T^. It is known 
that as the decay rate of hopping changes, one comes to a point at which there appears 
an infinite number of resonances, occurring at all length scales. The interaction of these 
resonances delocalizes the quantum states of the system. However, right at the transition 
point the density of resonances turns out to be low (provided that the coupling is weak), 
and one can construct a theory of the transition by treating the resonances as a kind of 
dilute gas. This makes it possible to construct a renormalization group theory (RG). We 
discuss the basics of the RG formalism in Sec. fV|, identify the quantities characterizing the 
strength of coupling, and then in Sec. derive an RG flow of the coupling parameters. 

The RG equation derived in Sec. has some interesting features. It has properties 
very similar to that of the Boltzmann equation, and can be studied by tools borrowed from 
statistical mechanics. In Sec. |VI| we flnd integrals of motion of the RG flow, derive an 
if— theorem, and obtain stationary distributions characterizing universal RG limit. 

The similarity between our RG equation and the Boltzmann equation has the following 
physical meaning. The RG procedure involves subsequent diagonalization of resonances 
ordered in energies, which is the same as the ordering in increasing spatial scales, i.e., in 
the RG time ^ = InR. Because of weak coupling, the only important kind of resonances 
are pair resonances which are well separated in the RG time domain. This situation is 
analogous to the statistical mechanics of a gas, where because of low density only pair 
collisions are important. Because of long collision time, one can neglect the correlations 
of subsequent collisions. Thus, everything about gas dynamics can be extracted from the 
mechanics of a two-atom collision, which is the essence of the Boltzmann equation. Similarly, 
in our problem the role of elementary collision events is played by pair resonances. The RG 
equation accounts for the change of the coupling parameters resulting from one resonance, 
and treats the resonances as uncorrelated events occurring randomly in the RG time. This 
analogy is quite fruitful: one can deflne an "entropy" obeying an if— theorem, and show 
that the RG flow leads to some analog of statistical equilibrium, in which the distribution 
of couplings takes a universal limiting form (which is the analog of Maxwell distribution in 
a gas). 

In Sec. |V11| we use the RG to study scahng properties of the states. For that, the RG 
is extended to include participation ratios of the states. We write down an RG equation 
for the flow of the participation ratio distribution, and by solving it determine the scaling 
exponent. The exponent turns out to scale with the coupling strength, Ze^/hv, and thus is 
small at weak coupling. This exponent has a geometric meaning of fractal dimension of the 
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space region occupied by an eigenstate. 

In Sec. [Villi we study transport in the momentum space, which in terms of Bloch states 
means angular diffusion. As it should be at the localization transition, the diffusion is of 
anomalous character, t ~ T>~^r'^. We find that the "diffusion constant" T> scales with the 
coupling strength. The dimensionless conductivity estimated from the Einstein relation 
(J = T>n is mush less than unity. Scaling properties of the scale invariant density- density 
correlator are related to the scaling exponent for the participation ratios derived in Sec |V11| . 

Finally, let us mention that the renormalization method employed in this work is not 
completely new. Similar ideas were used in the studies of one- dimensional quasiperiodic 
Schrodinger operators |]10|, and of localization in quasicrystals [llf. 



II. MODEL 

To be specific, we will focus on the problem of a particle moving in a 3— dimensional 
periodic crystal of scatterers: 

V{r) = Y^Uir-a), (2) 

agL 

where U{r) is a potential of one scatterer, and L is a 3— dimensional lattice with periods oi, 
02, and as: 

Oinii + a2n2j + a^usk , where Ua & Z . (3) 

The interaction U{r) has a power law singularity at small r, U{r) ~ \r\^'^, and falls rapidly 
at distances much larger than the lattice spacings Oq,, so that the sum is well defined. 
For the model we take 

Uir) = A e-IH (4) 

where the "effective radius" of the interaction is chosen to be of the order of the distance 
between the scatterers: ~ 1. 

We are interested in the spectrum and eigenstates of the Schrodinger operator 

+ V{r) (5) 



in the semiclassical range of energies: 



(6) 



Eigenstates in the periodic potential are Bloch states ipp{r) characterized by quasimomentum 
p in the Brillouin zone of the lattice L. 

Physically, the model (0)-(|D corresponds to the problem of the scattering of a fast 
particle by a crystal. Most interesting is the case of charged particles, like electrons or 
muons, because it turns out that the Coulomb potential singularity, f/(r) ~ - lies on a 
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critical line of a transition from regular to chaotic dynamics. However, it is tutorial to 
consider the problem with an arbitrary a. 

The problem is closely related to so-called quantum billiards. In the case of 

billiards, typically, one considers a particle moving within a rectangular, rhombic, or square 
billiard with a scattering potential inside and a hard wall boundary condition on the walls. 
In such billiards all nontrivial dynamics is due to the scattering by the potential, since the 
dynamics in an empty billiard is integrable. The problem of a hard core disk potential 
is known as a Sinai billiard. It is natural to generalize it to potentials with a power law 
singularity (|). 

The only difference between the billiard problem and our problem is in boundary con- 
ditions: hard wall for billiards vs. periodic in our case. It will be clear below that in both 
cases all results are the same. In fact, in our discussion of the mapping to a localization 
problem, as well as in the following treatment of localization, one can everywhere replace 
Bloch waves by standing waves without affecting any of the conclusions. Clearly, this is 
consistent with the mathematical equivalence of the billiard and the crystal problems that 
one has in the semi classical limit (^). 



III. MAPPING TO A LOCALIZATION PROBLEM 



Let us write the problem in the basis of Bloch plane waves e*'-''^"^^^^', where p is the 
quasimomentum in the Brillouin zone, and g is the vector of the lattice L', reciprocal to L: 

27r 27r . 27r , , 
g = — rill H n2] H nak , where no G z . (7) 

(1\ (I2 CL^ 

In this basis Schrodinger equation takes the form 

Ecg = E^c^ + J2 f^g-g'Cg'. (8) 
gYg 

where 

and Ug are the potential Fourier components = V^^ J e^^'^U{r-)d?r, where V = 010203 is 
the volume of the lattice unit cell. For the interaction M) one has 



^g ~ I^I.S-rrT/ / o, n . ,x(2-<t)/2 ' ) 



AttX T{2-a) sin((2-a)tan-ig//t 
|gF^ {K^g^ + 1 

which can be expanded at small and large g, correspondingly: 



where A„ = 47rAr(2 - a) sin(7r(T/2), = AnX r(3 - a). 

It is natural to think of the problem (P) as an Anderson localization problem defined on 
the sites of the lattice L', where E^s are on-site energies, and f/g-g' are hopping amplitudes. 
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Such an analogy is meaningful because the hopping is effectively local, since Ug falls at large 
|g|. In the problem we will find a delocalization transition at cr = 1 and any A, no matter 
how weak. 

Let us recall that in the conventional Anderson model with a hopping which is weak 
and short range, all states are localized. In our problem this is not so, because the hopping 
is actually long range: f/g ~ |g|^*-^^'^^ at |g| ^ k. In such a problem, even weak hopping 
delocalizes the states, provided the exponent in the hopping power law is less than the space 
dimension (see Sec. |I^ and Refs. [|12|,|13[). We will argue that the delocalization occurs over 
a shell of sites that have almost equal energies. In the momentum space, such sites are close 
to a "Fermi sphere" which has dimension two (see Fig. ^). Thus, at small A we have two 
phases: 

(i) 3 — 0" > 2, localized states; 

(ii) 3 — (T < 2, delocalized states. 

Localization in momentum space {a < 1) means that the wave function in real space is 
a linear combination of a small number of plane waves, or harmonics, i.e., it describes 
a quasiperiodic motion. On the other hand, the states delocalized in momentum space 
(cr > 1) correspond to a sum of many plane waves with almost equal wavelengths, but 
nearly isotropic distribution of wavevectors. In this case the dynamics is chaotic. 

To see why effective space dimension equals two, let us replace ?7g-g' in (||) by an infinitely 
long range hopping with a constant amplitude 

W = max{f/g} = B^K-^'^-'^^V-^ . (12) 

allg 

This problem corresponds to scattering by a 5— function potential, and to chaotic dynamics. 

Upon making the hopping infinite range all hopping amplitudes increase, and therefore 
the number of sites g over which the states delocalize should also increase. However, for the 
new problem the states can be easily written explicitly: 

where C'^'^^ is defined as the sum of all Cg's, 

C« = ^c», (14) 
g 

and the eigenenergy E^''-'> is found from the equation 

W 



Et;7I^ = 1- (15) 



g 



From that, one finds that the distribution of momenta in each state is isotropic, and that 
the sites g over which the z— th state will spread are all contained within the shell 

^(0 -^W <Eg< + -fW , (16) 

where 7 is of the order of one. Clearly, the quantity W just sets an upper bound on the 
delocalization shell in the problem (|^). 
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We are interested in the situation when W E, and hence the shell is very thin 
compared to its radius. Therefore, we can think of the shell as of an effective two dimensional 
"physical space" in which our localization problem is defined. The density of the on-site 
energies Eg can be calculated from the sphere \p\ = Po = V2mE surface area, 47rpQ, and 
the free particle density of states i^{E) = mpoV/2n'^h^ . The density of the on-site energies, 
taken per unit area of the shell and per unit energy, is 

n = VmpQ\27rh)-^ . (17) 

Another similarity with the localization problem is that the energies Eg can be treated 
as quasimndom numbers. To see this, one notes that EgS are given, up to a constant, 
by squares of the lattice vectors lengths. However, in a generic lattice L all vectors have 
different lengths, except for those related by a sign reversal g — g. This implies different 
distances to the constant energy sphere \p\ = po, and thus no occasional coincidence among 
E'g's. Thus, for now we assume the energies uncorrelated, and discuss consequences of their 
correlation later. 

To summarize, we argued that the problem (^-(^), from the point of view of the theory 
of localization, is equivalent to the problem 

Ecg = EgCg + ^g-g'Cg' • (18) 
gVg 

In the problem ([T8|) the sites g are uniformly distributed over a plane, and have uniform 
density n = VmpQ^ {27ih)~^ . It is assumed that in the semiclassical limit the number of 
lattice points around the constant energy shell becomes so large that one can replace the 
finite size sphere by an infinite plane without affecting localization properties. The nature of 
this approximation is similar to the one made in the conventional localization theory when 
one goes from a finite system of size much bigger than the localization radius to an infinite 
system. 

The energies Eg in (|T8|) are quasirandom numbers having unit spectral density, and the 
hopping amplitudes f/g-g' are defined by 

^7g-g' = A„e.|g-gr", (19) 

where 

\new = ATTXoidV'^r{2 - a) sm{na/2) and a = 3-a . (20) 

For convenience of notation, below we ignore the rescaling factor relating Xnew and Xoid, and 
denote Xnew by A. This will not cause any ambiguity because from now Xoid will not appear. 

Also, in the calculation done below we assume Coulomb potential for the scatterers, i.e., 
everywhere set o" = 1 and a = 2. The generalization to other potentials will be discussed at 
the end of the paper, in Sec. ^. 

Finally, let us mention that it is possible to generalize our problem to other space di- 
mensions. In fact, all arguments on the localization vs. delocalization behavior given above 
remain valid if one replaces the space dimension 3 by any d. In this case, the constant en- 
ergy shell in momentum space, in which the localization problem is defined, has dimension 
d — 1. The Fourier transform of a power law potential U{r) = Xr~" goes as Ug = X\g\'^~'^. 
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Then, by the locahzation criterion which requires that the power law in the hopping spatial 
dependence equals the space dimension, we get that the transition occurs at cr satisfying the 
condition 

a = d-a = d-1 . (21) 

This means that the Coulomb potential (a = 1) has critical singularity in all space dimen- 
sions (obviously, with d = 1 excluded). 



IV. DELOCALIZATION ARGUMENT: DIVERGING NUMBER OF 

RESONANCES 

We take advantage of the mapping to the localization problem, and from now on, until 
Sec. 1X1 we deal with the localization problem with long range hopping. For brevity, we will 



use d to denote the space dimension in which localization takes place, which corresponds to 
d — 1 in the above sections. The hopping power law exponent is a (equal to d^id — cr) . Let us 
recall that the character of localization depends on the nature of the hopping. For a short- 
range hopping the transition occurs at a particular ratio of the hopping to the on-site random 
potential. In contrast, for sufficiently long-range hopping, all the states are delocalized, no 
matter how weak the hopping is |^. More precisely, in a (i— dimensional systems with the 



|g|^" asymptotics of hopping localization can exist only for a > d (see |]T2[)- H en < d, 
then localization is destroyed, even at weak hopping. Thus, a localization-delocalization 
transition occurs at a = d. 

As a side remark, many real systems of interest do fall on the critical line a = d = 3 
(e.g., see [0). Among such systems are localized optical phonons in disordered dielectric 
materials coupled by electric dipole forces, two-level systems in glasses coupled by an 
elastic interaction, magnetic impurities in metals coupled by the RKKY interaction, etc. 

Let us briefly recall here the qualitative argument of Ref. showing that the states of 
the problem (|8|)- (p!oD cannot be localized. 

Consider two localized states located at g and g' having energies and E^'. They are 
resonance if 

l^^-g'l > \Eg - -Eg/ 1 (22) 



If the resonance condition (p^ holds, then eigenstates of the problem 



Ecg/ = Eg/Cg/ + f/g'-gCg (23) 

are not localized on g or g', but are essentially non-zero at both sites, g and g'. 

Now let us estimate the number of resonances and show that it diverges aX a = d. Let 
us takes one site g, and consider all sites g' which fall in resonance with the site g. The 
mean number of such sites is given by 



/ \gr-'d\g'\P[g,g'] , (24) 
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where P[g, g'] is the probabihty that the condition (p2D holds, and Hd is the surface area of 
a d— dimensional unit sphere. To estimate P[g,g'], we note that since the on-site energies 
Eg are uniformly distributed with the density n, one simply has 

P[g,g'] = |f/g-g.|n. (25) 

Substituting the probability ( pS]) with U^^gi ~ A|g — g'|~" into Eq. (|^), one gets an estimate 
for the average number of resonances: 

XuKd J |gr-""^rf|g'| . (26) 

This expression converges at large |g'| when a > d, and diverges when a < d. The divergence 
means that as one looks further away from a given site, there are more and more resonances. 
Obviously, in this case all the states must be delocalized. 

It is peculiar that the localization criterion does not involve the strength of hopping, but 
only the decay rate at large distances. In the standard localization problem, with a short 
range hopping, the localization transition can be reached by varying the hopping strength. 
Whereas in our problem, starting from the localized state, the transition can be reached 
either by increasing the hopping strength or by decreasing a. This opens up a possibility 
to study the localization transition at hopping amplitudes much smaller than the on-site 
potential variance. In this regime, due to relative smallness of the hopping, one can develop 
an accurate theory of the transition (see Sees. 0, |VlD . 

To see more clearly what kind of simplification one has in the weak hopping limit, let us 
look again at Eq. (^) and estimate the number of resonances at critical a = d. One has 

l^dn J Ig'r^lg'l = XriKd ln(i?max/-Rmm) , (27) 

where -Rmax is the size of the system, and -Rmin is a "microscopic" scale. The number of 
resonances (|27|) increases logarithmically with the system size, which indicates delocalization. 
However, if the hopping strength A is small, the effect of delocalization will be weak, because 
the resonances will occur rarely in the log— space. Another way to express it is to note that 
Eq. ( p7D means that the resonances have a uniform density over the log— shells: 

2^i?^in < |g - g'l < 2^^+^^^;^ , A; = 1,2,3,... (28) 

This suggests treating the resonances "shell by shell," i. e., a renormalization group approach. 
Since the density of the resonances distribution over the shells is small in A one can construct 
an accurate and self-consistent RG scheme. 



V. RENORMALIZATION GROUP THEORY. 



In this section we review the RG theory of resonances |T^. The reason an RG treatment 
is required in this problem is that when all direct resonances discussed above are combined 
together and diagonalized, there appear subsequent resonances, basically of the same nature 
as primary resonances. Then, after the next generation of resonances is treated, there appear 
more resonances and this process must be repeated (as shown schematically in Fig. ^). 
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To derive an RG equation an approximation will be made in which only pair resonances 
are taken into account, while all higher order resonances (triplets, quadruplets, etc.) are 
ignored. To motivate this approximation, let us discuss one property of resonances that will 
be basic for our approach. 

Let two sites g and g' form a resonance. Consider another site g" which is in a resonance 
with any of these two. According to Eq. (^) we can estimate the probabilities 

^{l- ^ 2) - An « 1 (29) 




^ ^ < 2 ~ An< 1 (30) 

|g-g'l / 

(here 2 can be replaced by any other number of the order of one). In other words, if the 
three sites g, g', and g" fall in a resonance, then one side of the triangle they form (say, 
|g — g"|) is much shorter than the other two (|g — g"| and |g' — g"|). In accordance with 
Eq. (pTl), one has 



( Ig-g" j - (^")" » 1 (31) 

Therefore, one concludes that pair resonances are much more frequent than triplets: the 
estimate (^) shows that the probability of a triplet resonance is of order of A^. Similar 
arguments show that the probabilities of finding resonances of k oscillators {k = 4, 5, 6, . . .) 
can be estimated as A'^~^. This should be compared with the probability of a pair resonance 
which is of order of A (see (^9]) ). We see that pair resonances occur about A~^ times more 
frequently than triplets, about A~^ times more frequently than quadruplets, and so on. This 
justifies keeping only pairs in the RG. 

The RG procedure involves the following steps. We truncate the Igl^'^— interaction at 
some Rq, i.e., put f/g-g' = for all pairs (g, g') such that |g — g'| > Rq. Let us assume 
we know exact eigenstates for the truncated hamiltonian (call them Rq— states). Then, to 
eliminate next renormalization shell, we replace Rq by Ri such that 

R 

i?i > i?o , but An logs < 1 (32) 

ito 

Let us take the -Ri— states and represent them as linear combinations of the -Ro~states. 
According to the above discussion, all of states are either single i?o— states or resonance 
pairs of -Ro~states (we neglect by triple and all higher order resonances). 

The consistency of the RG procedure relies on the fact that the localization radius of the 
i?— states is less than R. One notes that by going from i?o— states to states, at weak 
hopping, the localization radius cannot exceed Ri. In other words, on the scale R the states 
will look like a superposition of several states localized on a scale much smaller than R. 

Practically, this means that at all R one can work in the basis of localized states. In 
the RG scheme, when considering the -Ri— states as a result of interaction (resonance) of 
the i?o— states we assume that the resonating i?o— states are far apart compared to their 
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localization radii. This enables to keep track of the coupling between the states, ignoring 
the details of the spatial structure of wavef unctions. 

The RG flow preserves the form of the Hamiltonian, but modifies the hopping amplitudes. 
We will see below that the change of the amplitude ?7g-g' of hopping between g and g' is 
given by a multiplier that factors into Og/Og. The quantities Og have a meaning of coupling 
parameters generated by the RG flow. Mathematically, at each R the truncated Hamiltonian 
has the form 

Ecg = EgCg + Ug-g>agag' Cg/ , (33) 

\s'-s\<R 

where the quasirandom energies Eg and the hopping f/g-g' = A/|g — g'|" are the same as in 
the original problem (^-(p!OD, and the new quantities are the parameters Og characterizing 
the strength of coupling between the states. Initially, all Og = 1. However, the RG procedure 
generates a non-trivial distribution of the parameters Og. 

To see this in more detail, let us write down all necessary quantities for two states (at 
R = Rq) located at g and g'. They interact according to 

-Z^Cg — -/^'g Cg I ^^g ' 

ECgi = EgiCgi Ug'—gCLgiCLg Cg , 

Two eigenstates c+ and c~ are defined by 

= cos ^ Cg + sin 9 e**^ Cg/ , c~ = — sin 6* Cg + cos 9 e**^ Cg' , 

with 

/ 

tan 6' = + 1 — r , r = ^ — — — - , (p = arg(agag/) . (36) 

2|(7g_g'CigCig' I 

The energies of the states (^) are 

E± = l [Eg + Eg, ± ^(Eg-Eg02 + 4|[/g_g,agag,p) (37) 

The transformation rule for the parameters Og can be found from the relation 

o^gCg + Og'Cg/ = a^c^ + a^c 

which means that 

a'*' = cos 9 ag + sin 9 e~^'^ Og/ , a~ = — sin 9 ag + cos 9 e~^'^ Og/ (39) 

The transformation rule (|39| ) for the coupling parameters Og shows that even when initially 
all flg = 1 the resonances lead to non-unit values of Og, which justifies including Og's in the 
RG. 

As a matter of fact, the parameters Og are the only important ingredient of the RG. Other 
characteristics of the states, such as their energies, turn out to be irrelevant. The reason is 
that when a resonance of two states forms (say, at the scale R) the energies of the states split 
by something like |?7g-g'| ~ -R~". (Since at the scale R the resonances are typically spaced by 

12 



(34) 



(35) 



(38) 



|g_g'| ji '^ This splitting ensures that at larger scales these states never fall in resonance. 
Because of that, the energies Eg can be taken as uncorrelated quasirandom numbers with 
uniform distribution, and their change under RG transformations can be ignored. 

Now, let us proceed with deriving an RG equation for the distribution of Og's. In doing 
it we emphasize the analogy with the theory of Boltzmann kinetic equation. Boltzmann 
equation for gases is usually derived assuming the absence of correlations of subsequent 
collision processes, resulting from large mean free paths of molecules in a gas. Besides 
providing grounds for a probabilistic approach, the largeness of mean free path enables one 
not to take into account triple and other multiple collisions. Similarly, in our RG dynamics 
there is no repetition of resonances (collisions) of the same states, because of the energy 
splitting. Also, small density of the resonances in the log— spa.ce, analogous to large mean 
free path in gases, makes triple and higher order resonances unimportant. 

As a first step of deriving the RG equation, let us consider how the distribution of a^s 
changes by going from i?o— states to states. The distribution /(a) is defined by the 
probability density 

dP = f{a)dadd . (40) 

Under the change Rq Ri the distribution /(a) changes to f{a). The change is due to 
resonances ("collisions") occurring during the RG time interval lni?o < C, < lni?i. The 
difference f (a) — f (a) can be obtained by integrating 
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6 (a — a^^ + 6 (a — a ^ — 6 {a — ai) — 6 {a — 02) (41) 



over 

1) f (ai) daiddi] 
3) f (^2) da2da2; 

3) nd'^g , where g = gi — g2 and n is the concentration of sites given by Eq. (|l] 

4) dE , where E = Ei - E2. 

Here the subscripts "1", "2" label the states in a resonance pair. It is convenient to 
introduce a new variable r instead of E, defined by the second equation in (PB[). The 
rotation angle is related with r by tan 6 = y/r"^ + 1 — t. The usefulness of the variable r 
becomes clear from the identity 



d^gdE = 27r|g|(i|g| 2\aia2Ug^-g^\ dr = 47r|g|^|aia2f/gi-g2 1 |g| '^d\g\dT 
= (47r|gi - g2\^\aia2Ug^^g^\j d{\n\g\)dT 



(42) 



where the product |gi — g2p|'5ia'2f^gi-g2l does not depend on |gi — g2|. (Because t^gi-ga = 
A|gi— g2|~^.) Thus, we perform an integration over gi— g2 in the domain Rq < |gi— g2| < -Ri. 
The result is 

f {0) — f {0) = ln(i?i/i?o) j dTddidaidd2da2 \a2\ f (oi) / (02) 

X 6 (^a — a~^^ + 6 (^a — a'^ — 6 {a — ai) — 6 {a — 02) (43) 

According to the above discussion of the probability of pair resonances, the RHS in ( ^3]) is 
of order of An In (Ri/Rq), which is a small number (see (|32|) ). Then it is standard to replace 
Eq. ( ^3|) by a differential equation with respect to the RG time ^ = In (R): 
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= Xn J dTdaidaida2da2 \ai\ |ct2| 
X / (ai) / (02) 6 (^a — a^^ + 6 (^a 



6 (a 



5 [a 



(44) 



This equation is the main result of this section. Our task now will be to look for solutions 
/(a,^) to the problem ( ^l] ) with the 5— function initial condition /(a,^ = 0) = 5{a — 1) 
corresponding to the "microscopic" distribution of the parameters a,. 



VI. SOLVING THE RG EQUATION 

The similarity with the Boltzmann equation makes the dynamics ( ^ ) very simple. The 
problem ( ^If ) has integrals of motion (analogous to the conservation of energy and momentum 
for Boltzmann equation), and also it obeys an if— theorem. Relying on these properties we 
analyze the asymptotical behavior of / (a, C) at ^ — > 00. 

First, let us study the integrals of the dynamics (H). The simplest integral is 



(|a|^) = / dada\a\^ f{a,$) 



Let us prove that 
d 



(|ar) = o 



or 



j |a|^/(a, ^) dada = J |a|^/(a,0) dada 



Proof: 



d_ 



(|ap) = n dr daidaida2da2f {ai) f {a2) 



+ 



- ml 



r2 



(45) 



(46) 



(47) 



But 





2 




2 


a+ 


+ 


a 





cos 6 tti + sin 9 e ^'^ a2 



+ 



sin 9 tti + cos 9 e a2 



|ai|^ + |a2|^ (48) 



(unitary transformations preserve the norm). Consequently, the RHS of Eq. (^) vanishes. 
QED 

Besides there exist other invariants of Eq. (^4]). All averages of the real and 

imaginary parts of a = a' + ia", 



{a'a") 



(49) 



are conserved by the same argument of orthogonal transformations preserving the norm. 

An interesting issue to discuss is whether the conservation of the quantities (^) is an 
exact or an approximate result. One might worry that the conservation fails when higher 
order resonances are taken into account, in addition to the interacting pairs. However, the 
conservation laws remarkably survive. To see why that is, let us consider a resonance of k 
states. The Hamiltonian is diagonalized by an orthogonal k x k matrix Ua/3, (a,/3 = 
1, . . . ,k). One notes that the parameters Oq, are transformed by the the matrix U~^, and 



then the conservation of the quantities (|49|) follows from orthogonality of U 
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Another property of the problem is an analog of the if— theorem. Let us define 
entropy of the distribution f{a,^) as 

[f] = -J In (|a|/ (a, 0) / («, dada (50) 
The if— theorem states that [/] grows monotoneously as function of the RG time ^: 

[/] > . (51) 
Proof: Using the identity = ^1 = ^/ / (a, ^) ddda = J ^/ (a, ^) (ia(ia one finds 



laaa 

|«i|/(«i) |«2|/(a2) 



dH f d f 

— = - J In (|a|/) —ddda 



n dr 



J J daidaida2da2\ai\f (ai) \a2\f {a2)\n. 



\a+\f{a+)\a-\f{a') 



(52) 



Another useful identity 



J dr J J \ai\ f (ai) \a2\ f {a2) daidaida2da2 ~ J J l*^^!/ ('j^^) k 1/ ('^ ^da^da^da da 



= 
(53) 



follows from the fact that the variables ai, 02 and a+, a are related by an orthogonal 
transformation. Combining the identity ( ^3|) with (|52D one gets 



^ = n J dr J J daidaida2da2^ (a+) ) (-^ln(X) - X + 1) , (54) 

where X = ($ (oi) $ (a2))/($ (a+) $ (a")), and $ {u) = \u\f (u). 

Finally, the quantity X ln(X) — X + 1 is non-negative. From that, the if— theorem (|5TD 
follows momentarilly. QED 

Note that although our entropy differs from Boltzmann entropy by (In |a|), it does not 
affect the proof of the ii— theorem which follows the standard derivation (cf. Ref. [|T5[]). 

The problem (^) has stationary distributions given by 

fcia) = j^exp{-aaGapai3) , (55) 

where ai = a', a2 = a", and G is a positively defined symmetric 2x2 matrix. 
Proof: If a solution f{a,^) of Eq. ( PI) does not depend on ^ then ^ = for it. Since 
Xln(X) — X + 1 = vanishes only for X = 1, according to Eq. (Q) for such a solution 
$ (oi) $ (02) = $ (a^) $ (a~) for all 6 (hence, for all r). This is possible provided that 

/(a)=/G(a) • (56) 



Given Ga/3, the normalization constant A in (|55|) should be determined from the condition 
/ / (a) ddda = 1. <5EI) 
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Stationary solutions fa (a) are analogous to the Maxwell distribution which appears 
in statistical mechanics as a stationary solution of Boltzmann equation. An important 
distinction is that in the Maxwell distribution there is only one free parameter (temperature) 
while the distributions /g (a) are characterized by 3 parameters Ga/3 (a,/? = 1,2; a < 
/?). Like in the Boltzmann equation theory, the number of parameters in the stationary 
distribution is equal to the number of conserved quantities (^91). 

Continuing the analogy with statistical mechanics, let us prove that our entropy is max- 
imal on the stationary distributions (|55|). More precisely, let us take all functions / (a) such 
that 



/ (a) dada = 1 , J aadfsf (a) dada = Gap , («,/? = 1, 2) 



(57) 



where ai = a', 02 = a", a = a' + ia" . Then always H[f] < H [fa]-, where fa {a) is defined 
by (HD; [/] = H [fc] only if / (a) = fa (a). 

Proof: Let us consider the first and the second variations of H[f\ constrained by ( ^7]) and 
by the normalization condition: 



H [f] — \ j f (a) dada — ^ /i^^ J aadpf (a) dada 

Sf (a) dada 



-ln(|a|/ (a)) - 1 - A - J2 l^c^paaap 

a/3 



(58) 



H [f] — \ f (a) dada — ^ /Xa/3 / cta^'pf {0,) dada = — — TT^ — dada < 0. (59) 

■J a/3 •'J W 

Since the second variation is negative, maximal value of H [/] corresponds to 



/ (a) = 777 exp I - (1 + A) - ^ lia^a^ap | (60) 

a/3 



a 



The relationship between A, /ia/3 and A, Gap is straightforward (compare (^) with (|^) ). 
QED 

The asymptotic properties of solutions of Eq. (^^ are thus rather simple: each solution 
converges to one of the stationary distributions fa (a)- The parameters Gap are completely 
determined by second moments of the initial distribution / (a,^ = 0): 

Gap = J CLaO-pf (o, ^ = 0) ddda (a, /5 = 1, 2, a < /?) . (61) 



To conclude, the RG dynamics ( ^41) has a 3— dimensional family of fix points param- 
eterized by symmetric positively defined 2x2 matrices. The dynamics in the space of 
distributions /(a) is shown scematically in Fig. ^. 

The parameters are determined by the initial distribution, which means that the coupling 
strength in this problem is not renormalized. Thus the critical behavior such as scaling 
of the states will depend on these parameters. This situation should be contrasted with 
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the ordinary localization transition at short-range hopping, where the critical behavior is 
believed to be completely universal, with no additional parameters. 

Let us mention that for the particular problem of Coulomb scatterer we are interested 
in the distribution /^(a) has the form 

f{a') = A\a'\-^exp(^-Aa'^) , (62) 

i.e., it does not depend on the imaginary part a". The reason is that in the initial distribution 
all parameters a are real: f{a)\^=o = 6^^\a-l). The non- generic form of matrix Gaf3 resulting 
from the fact that all a's are real leads to a complication: the stationary distribution (|62[) 
is not normalizable. However, since the divergence of / f{a')da' is only logarithmic (i.e., 
weak), the results in this case are the same as for generic Ga^. The only subtlety is that 
the constant A in ( |52[ ) is renormalized, at a rate proportional to In^, which leads to the 
distribution slowly flowing to more and more singular distributions. However, all results 
derive above still hold, including the absence of RG correcions to the coupling strength 
((a')^). Because of that, the conclusion about non-universal RG fix point holds in this case 
as well. 



VII. WAVEFUNCTION SCALING: RG FOR THE DISTRIBUTION OF 

PARTICIPATION RATIOS 

In this section we study renormalization of the eigenstates. The hierarchical structure 
of the states (resonances of resonances of resonances ...) emerging from the RG dynamics 
suggests that they have fractal properties. To study the scaling behavior of the states we 
will consider their participation ratios. For a normalized state c^*-* its participation ratio p'-*-* 
is defined by 

p« = Eic«r (63) 

g 

In the theory of localization participation ratios are used as a measure of the degree of local- 
ization. By studying scaling of participation ratios one can characterize fractal properties 
of the states. 

Below we study how participation ratios are changed by the real space RG of the states 
described in Sec. For this purpose we need participation ratios p^^ of the i?— states, which 
have been defined as true eigenstates of the system with the interaction truncated at the 
scale R. To incorporate the participation ratios in the RG, one can rewrite the RG formalism 
for the distribution f{a,p), including both coupling strengths and participation ratios. In 
fact, the i?— states are characterized not only by a^^^ and p^'^\ but also by eigenenergies 
and positions. However, as we already discussed, the correlations between the energies and 
positions are not important, and thus at each step of the RG one can assume uniform 
uncorrelated distribution over energy and in space. 

As it was argued above, the RG fiow in the leading order is determined by pair resonances. 
Let us discuss how pair resonances affect the participation ratios, and derive an RG flow 
for the distribution f{a,p). For that we consider the resonance pair (^) and flnd for it 
the participation ratios of the states "+" and "— " as functions of Pi and pj. Since the 
positions of the states i and j are well separated, one simply has: 
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cos^ 9 Pi + sin^ 6 pj, 



p 



sin^ 6 Pi + cos'* 6 pj . 



(64) 



This change of the participation ratios must be considered together with the transformation 
of the parameters Oj, aj given by Eq. (|39D. Other steps of the derivation of the RG flow 
for / {a,p) are similar to those described in Sec. 0, and the resulting RG equation is the 
following: 



d_ 



f{a,p) = y dTdaidaida2da2dpidp2 \ai\ \a2\ / (oi, Pi) / (025^2) 

X [5 {a — a^^ S (p — P^) + 6 (^a — a^^ ^ (p — 
- 5{a-ai)5{p- pi) -5{a-a2)5{p- P2) ] • 



(65) 



Of course, from the flow equation (|65|) one can return to Eq. ( ^41) by using the relation 
between the functions / (a) and / (a,p): 



f{a,p)dp 



(66) 



By using the RG flow for the distribution f{a,p) given by Eq. ( pSD one can study the 
.^—dependence of various moments of p and of a, such as (p^) or as well as correlators 

{p^\a\"^). Here we concentrate only on the distribution of the participation ratios, assuming 
the distribution of the parameters a to be stationary. In other words, we will be interested 
only in the scaling limit of the problem. 

Let us focus on calculating the average participation ratio 



(P) = J J J dadadp pf {a,p) \^=inB^ ■ 
To facilitate the analysis we introduce a new function s (a) defined by 

1 

pf {a,p)dp . 



s [a] 



(67) 



(68) 



over p: 



The RG flow for s{a) can be obtained by integrating Eq, 

d f _ _ 

— s (a) = ^ J dTdaidaida2da2dpidp2 |a2| / (cn,Pi) / ((^2,^2) 
X [p'^6 {a — a^j + p^6 {a — a^^ — pi6 {a — ai) — P26 {a — 02) ] . 

The integrals over pi and p2 in Eq. (^) can be calculated by using 
d 



and (16 



—s (a) = n I dTdaidaida2da2 \ai\ \a2\ s (oi) / (02) + s (02) / {ai] 



where 



and 



= cos^ 6 5 (^a — a"*"^ + sin^ 6 6 (^a ~ a ^ — 5 (a — ai) 



(69) 



(70) 



(71) 



sin 6 6 ia — a 



cos 6 6 (a — a 



6 {a — 02) 



(72) 
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In the scaling limit ^ oo, as we demonstrated above, the distribution / (a) converges 
at large ^ to a stationary distribution fcio)- Since we are interested in the asymptotical 
behavior of s (a) at oo), the function / (a) in Eq. (^) can be set equal to the stationary 
distribution fa (a) given by Eq. (|55|) . After that we get a linear integral equation for s (a) 
with a kernel idependent of S,. 

With that in mind, one can view Eq. ( [70D as 



^^s{a) = -B{s{a)) , (73) 

where i? is a constant (^—independent) linear operator which acts in the space of functions 
s (a) . Asymptotical behavior of s (a) at ^ — >■ oo is then given by 

s (a, - exp (y-ni) So (a) , (74) 

where fi is the lowest eigenvalue of B and Sq (a) is the corresponding eigenvector: B (sq (a)) = 
fj,So (a). Then, according to Eq. (j^, the scaling behavior of (p) is given by 

(p) = J dadas (a) \^=inR ^ -R"'' ■ (75) 

Therefore, the scaling exponent of (p) is related to fi, the lowest eigenvalue of B. 

In principle, if the stationary distribution fc (a) is known, fi can be determined (per- 
haps, numerically). However, here we are interested in estimating /x only by the order of 
magnitude. 

To estimate fi, let us note that the kernel in Eq. (|7^) is of the order of 

n{\aia2\) — n{a^) = Xn . (76) 

Hence, given n and A, 

/i = 7An , (77) 

where 7 is a constant of the order of one. 

It is tutorial to derive the result ([77| ) by another method. Using Eq. (^) one can write 



-/i = — In {{p)) l^^oo = {p) ^ J -Q^s (a) dada 

= n{p)^^ dTdaidaida2da2 \ai\ \a2\ s (ai) /g (02) + s (02) /g (o-i) [••]2 > (78) 



where 

[..]^ = [..]^ = sm^e + cos^e-l. (79) 



But sin 9 + cos 6 — 1 = —2 sin 6 cos 6 < 0, and hence Eq. ( |78D can be rewritten as: 

— jj = —An J dr sin^ 9 cos^ 6 J dddafda) (^J ddda\a\s{a) ^ (^J dddas{a)^^ (80) 
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By estimating the RHS of Eq. ([80| ) by the order of magnitude, one again obtains the result 
([77|) and in addition proves that /i is positive. 

The exponent fi has a simple geometric meaning. For a typical i?— state, it gives the 
fractal dimension of a region occupied by this state. Hence, we find that in this problem the 
fractal dimension of the states depends on the coupling strength An. 

The argument presented above for the mean participation ratio, after proper modifi- 
cation, is also valid for any moment of the participation ratios. Hence, all other fractal 
dimensions must scale with An the same way as (p). 

In the localization problem it is sometimes of interest to consider the distribution of all 
moments of the wavefunction: 

p{m)=J2\cfr. (81) 
g 

The participation ratio studied above corresponds to m = 4. The distribution of p{m) can 
be used to study fractal dimensions of the states. 

It is straightforward to generalize our RG formalism for the quantities p{m). For that, 
one simply notes that for the pair resonance the change of p{m)^s is written as 

p~^{m) = I cos^l™ Pi{m) + I sin 6'!™' pj{m), p~{m) = \ sin^^P Pi{m) + | cos^P Pji"^) ■ 

(82) 

Putting it together with the transformation ( P^ leads to the RG flow for the distribution 
f{p{m),a) of the same form as Eq. (|65D, but with p^ given by Eq. (|82D. The analysis of 
this RG equation is similar to that presented above, and so are the results. The scaling 
exponents for the quantities p{m) are of the order An. This confirms our conclusion about 
fractal dimension of the space region occupied by an eigenstate being of the order of An. 



VIII. TRANSPORT: ANOMALOUS DIFFUSION AND SCALING 

In this section we study time-dependent transport. Particularly we are interested in the 
density correlation function 

K{g- g', t) = Y: (e-^(^-^^)* c«4''^4' Cg' ) , (83) 

where c^*^ are eigenstates, Ei are their energies, and (...) stands for averaging over the 
on-site disorder. 

We will argue that the dynamics is critical (anomalous diffusion) and write down a scale 
invariant expression (|89D for the correlation function K{g,g',t). Then, by using the RG for 
the participation ratios presented in Sec. |V11| we will derive scaling properties of the function 



K. 

Let us begin with a qualitative picture of propagation of an excitation in the system. 
Let the particle start at t = from one of the sites (say, from the site 1 in Fig. ^). After 
some time T it hops to its nearest resonance neighbor, site 2. The hopping time T is of the 
order of inverse resonance splitting, which, according to Sec. |rV|, can be estimated as 
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T~ |^+-^_|-^ ^ Iri-ralVA . (84) 

After the time of the order of T the particle is locahzed not only at the site ri, but both at 
ri and at r2, since both modes c"*" and c~ are excited. 

Later, after some larger time T' ^ T new site will be involved in transport. Each of the 
two excited states sets in motion other states which themselves can be considered as a result 
of interaction of several states during the time interval < T'. This is shown in Fig. where 
new resonance states participating in the transport are linear combinations of C3, C4 and of 
C5, Cg, respectively. Obviously, L' = |ri — rQ,| |r2 — rd ^ |ri — r2|, where a = 3,4,5,6. 
Thus, at times t which are longer than T' but shorter than the time necessary for exciting 
new states, the excitation is mainly carried by four states. According to the estimates made 



in Sec. |^ one has 

L'/L- In-r^l/ln-ral ~exp(const/(An)) , T'/T ^{L' /if, (85) 
where a = 3, 4, 5, 6. 

Even later, after some time T" (T" ^ T') new states are involved in transport. The 
corresponding sites are separated from the initial site by a large distance L" such that 

L"/L' exp (const.A-^) , T"/T' ~ {L"/Lf . (86) 

At this stage, the number of modes carrying the excitation again approximately doubles. 

Let us remark here that the excitation evolution described above should be understood in 
a statistical sense. For a particular state, there is no exact length scale at which the number 



of resonances doubles. As we saw in Sec. IV, the hierarchy of resonances is characterized 



by density uniform in the log of spatial scale. The same is true for the time hierarchy of 
resonances. 

To summarize, the larger the time, the more states are excited. Asymptotically, after a 
large time T about 

2cAnlog{T) rpcXn (8T) 

states are involved in the transport (here c is a number of the order of one to be determined 
later) . 

There are two conclusions one can draw. First, the RG picture of the states hierarchy 
in space is in a one-to-one correspondence with the hierarchy in the time domain. The 
sequence of times at which new states are involved in dynamics is determined by the spatial 
hierarchy of resonances. 

Second, according to the above discussion, after a long time T the excitation spreads over 
a region of size L (T) ~ T^/^. This means that the dynamics is (asymptotically) invariant 
under rescaling 

{r,t)^(Zr,ZH) , (88) 



where the rescaling factor Z is an arbitrary number. 

In particular, the correlation function K (r, t) must be invariant under the rescaling 
Thus we have 
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at large t, |r|. The function -F(x) is normalized, 

J F (x) ci^x = 1 , 



(89) 



(90) 



due to probability conservation. The particular form of -F(x) depends in some universal way 
on the distribution of couplings /(a). 

Let us note here that the simple scaling relation (^) is a direct consequence of the 
existence of stationary fix points fcio,) of the RG flow for /(a) given by Eq. (|4^) . Indeed, 
the coupling of complex modes composed of many single-site states is determined by the 
distribution of the parameters Oi given by the limiting distribution fc (a) = lim^_^oo / (a, 0- 
Because this limit exists, the coupling of complex i?— states asymptotically scales as R~^. 
The time of hopping scales as inverse coupling, which leads to the space-time scaling relation 

There is a relation between the behavior of the scaling function -F(x) and small x and 
the participation ration scaling exponent /i discussed in Sec. [VIII| (see Eqs. ([75|) , ([77D). It 
turns out that 



|-(2-m) 



at 



Ixl < 1 . 



(91) 



To prove this relation, let us consider a finite system of size R and take the density correlator 
K given by ([83| ) at times t < R^. By the order of magnitude, K can be estimated by removing 
from Eq. (^) all terms with non-equal energies Ei ^ Ej. (This is correct because typically 
\Ei — Ej\ is larger than level spacing, which is of the order R~^, and hence the phase factors 



{Ei — Ej)t in (|83D are large, which leads to random relative signs of all terms with i ^ j, 
and leads to mutual cancellations.) 

On the other hand, after the terms with i ^ j are removed, one gets the following 
estimate: 



Ei4"rkgr 



(92) 



Now, let us set 



and sum over all 



Ei^(g,g', t^i?^) = E 



g =g 



where p*^*-* are the participation ratios. The right hand side of Eq. 
(see Eqs. ([75D, (fTTD). This is consistent with the scaling form 



K{g,g',t)= t-'F{\g~g'\/t 



(93) 

I) is of the order R^~^ 
(94) 



only if the function F satisfies Eq. (0). 

Also, let us remark that in an arbitrary space dimension d the critical law of hopping, 
f^g-g' ~ |g ~~ g'l~*^) leads to anomalous diffusion: T ~ R'^. The reason is that, as it is clear 
from the above discussion, the scaling transformation (BBI) changes to 
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(r,t)^(Zr,Z'^t) , 



(95) 



and the scale-invariant function -F(x) is now defined by 

K{r,t) =t-^F (t-^/'^r) . (96) 

(/ F (x) (i'^x =1). The anomalous diffusion follows directly from (p^): 

(|r|2) = J \r\^K (r, t) d'^r ^ t'''' j |xpF (x) d'^x . (97) 

In the problem of three-dimensional Coulomb scatterers, the space in which we have a 
localization problem is two-dimensional (constant energy surface in momentum space), and 
thus in this case there is no distinction between anomalous diffusion and ordinary diffusion. 

To compare the results with the standard localization theory assuming short-range hop- 
ping, let us point out that diffusion constant in owr problem is of the order A. Dimensionless 
conductivity can be found from the Einstein relation a = nD to be of the order An <^ 1. 
Thus in this problem one has critical behavior with a <^ 1. On the other hand, in the 
delocalized phase of the system with short-range hopping one always has a < 1. 



IX. SCATTERING OF CHARGED PARTICLES IN CRYSTALS 

Let us return to the problem of charged particle scattering in a crystal, and consider the 
dynamics of a fast charged partcile (electron or muon) in a real crystal. Since atoms cores 
are just 1/r— scatterers, one expects that the effects of weak chaos we discussed above will 
take place, and lead to delocalization in momentum space, i. e., to the absence of Bloch 
waves. Of course, from a practical view point, there is always some inelastic scattering, 
which complicates observing the effects we are interested in. However, in the following we 
ignore such processes, as well as all possible effects of crystal boundaries. 

Let us estimate the dimensionless coupling constant An. From ([TtD and (11) one has 



An = — , (98) 

where Z is the nucleus charge and v is the particle velocity. The coupling constant is of the 
order of one when the velocity v is of an atomic scale. 

From (^) it is clear that with electrons one never has a strong coupling situation for 
a semiclassical problem. The reason is that if an electron is moving with velocity of an 
atomic scale, v = Ze'^/h, or slower, the number of wavelengths in the region where the core 
potential is not screened is of the order of one. On the other hand, if the velocity increases 
so that the wavelength becomes small, one has a well defined, but not a very interesting 
semiclassical problem, because the coupling ([98| ) becomes very weak. Therefore, channeling 
of fast electrons (or positrons) is not a suitable situation. 

On the other hand, one gets into the desired regime much easier with heavier particles, 
of which the best studied is channeling of positive and negative muons. If a muon travels 
with velocity of an atomic scale Ze'^/h, so that the coupling (|98|) is of the order of one, its 
momentum is very large: 
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P/i — Patomic 



me as as 



(99) 



Correspondingly, the particle wavelength will be very small compared to the core potential 
radius as, and thus the problem is semiclassical and our results are applicable. 

The theory we described above predicts delocalization in momentum space. In terms 
waves scattering it means that a plane wave, after multiple scattering on the core potentials, 
will give rise to many other plane waves, and eventually there will be almost no memory 
preserved about initial plane wave. Following the discussion of Sec. |IV|, the number of plane 
waves which are in a resonance with given plane wave can be estimated as 



At large p^/patomic one can have many resonances between plane waves even for a weak 
coupling situation. An < 1. In this case, the theory presented in this paper will be accurate. 

There are several simple predictions which can be tested experimentally. It follows from 
our discussion of resonances that the number of plane waves into which initial plane wave is 
scattered, although large, is much smaller than the total number of plane waves with energies 
close to the initial energy. The origin of this effect is in the critical nature of delocalized 
states for Coulomb potential. The diffusion over constant energy surface in momentum space 
leads to a nontrivial angular correlation function, with a power law singularity (^1]) at small 



Our discussion of the localization in momentum space focused on the Coulomb scatterers 
which have marginal singularity and give rise to critical states in momentum space. It is easy 
to see how the RG results will change for a general power law singularity, U{r) = Xjr" . The 
hopping in momentum space, given by Fourier transform of U{r\ will be more short-range if 
0" < 1 and less short-range if cr > 1. Hence, at a > 1 the system is always in the delocalized 
phase, no matter what A. On the other hand, for cr < 1 the RG will scale coupling down to 
zero, and we can conclude that for sufficiently small A (i.e., when the RG is applicable) all 
states are localized. 

To represent this graphically, one can draw a schematic phase diagram in the plane (a. A) 
(see Fig. ^). The region marked by I corresponds to the standard delocalization transition 
for short-range hopping (large a). The region marked by II corresponds to the transition 
at small A studied in this article. One can note that the critical behavior in the cases I 
and II is quite different: all scaling relations in the case II depend on the dimensionless 
coupling strength which does not change under the RG flow, whereas in the case I there are 
no parameters. This means that one cannot go continuously from one scaling regime to the 
other, and hence we conjecture that there must be a point on the localization-delocalization 
phase boundary at which the critical behavior changes. 

Another reason to believe that there is such a point is that in the RG theory presented 
above the transition occurs at the same a = \ for any A (which must be sufficiently small). 
So the phase boundary in Fig. ^ has a part that goes straight up. This suggests that the 




(100) 



angles. 



X. GENERAL PICTURE 
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special point at which the critical behavior changes must simultaneously be the point C 
where the phase boundary turns away from the a = 1 straight line. 

One more reason supporting this phase diagram follows from recent paper by Mirlin and 
Fyodorov (see Ref. ll^). They studied a one-dimensional problem with power-law hopping 



in the regime of strong hopping where one can use the supersymmetric sigma-model. In 
the {a, A) plane this corresponds to large A. In this region, a delocalization transition was 
found at certain critical cr. It is natural to make a connection with our results by assuming 
that the point C is a tricritical point at which all phase boundaries match. 

To summarize, in this article we consider the problem of Bloch states for a fast particle 
moving in a crystal of scatterers with power law singularity. This problem is shown to 
be equivalent to a localization problem in momentum space which we study using an RG 
previously developed for the localization with long-range hopping. We find a localization 
transition with the critical singularity being that of Coulomb potential. In the localized 
phase the particle state is close to a single plane wave, whereas in the delocalized phase 
the state is a superposition of many plane waves with wavevectors spread uniformly over a 
constant energy shell in momentum space. In other words, the localization transition is the 
transition between integrable and ergodic dynamics. We believe that the relation between 
delocalization in momentum space and the transition to chaos is very generic, and will be 
useful beyond the problem of power-law scatterers. One can hypothesize that understanding 
of Quantum Chaos in general is in some way a problem of finding an equivalent in localization 
theory. 
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FIGURES 




• •••••• 

FIG. 1. A sphere of constant energy is shown in momentum space, together with a shell of 
the dual crystal lattice sites (0) which are close to the sphere. Delocalization of Bloch states in 
momentum space occurs within this shell. Effective dimension of the space in which we consider 
the localization problem is given by the sphere dimension, i.e., equals two. 




FIG. 2. Hierarchical structure of two eigenstates is shown: a) (((l+2)+(6+5))+((8+9)+ 7)) ; 
b) (((l+2)+(3+4))+((10+ll)+(12+13))) ; where we denote by (n + m) a resonance pair formed 
by the states m and n. The time hierarchy of excitation transport is given by the same diagram. 
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FIG. 3. Schematic view of the RG flow given by Eq. (|4^). The space of functions /(a) is 
represented by a bundle of invariant surfaces labeled by the integrals Gap = {aadp), {a,f3 = 1,2, 
a < f3). Restricted on each of the surfaces the system has one attracting fix point (given by 
Expr. dH)). 
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FIG. 4. Phase diagram in the (cr — A) plane. (Both a and A vary between and oo.) The 
phase boundaries are marked corresponding to the delocalization by long-range forces studied in 
this work (region IT) and to the conventional Anderson transition at short-range hopping (region 
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